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In graphene spintronics, interaction of localized magnetic moments with the electron spins paves 
a new way to explore the underlying spin relaxation mechanism. A self-assembled layer of organic 
cobalt-porphyrin (CoPP) molecules on graphene provides a desired platform for such studies via 
the magnetic moments of porphyrin-bound cobalt atoms. In this work a study of spin transport 
properties of graphene spin-valve devices functionalized with such CoPP molecules as a function 
of temperature via non-local spin-valve and Hanle spin precession measurements is reported. For 
the functionalized (molecular) devices, we observe a slight decrease in the spin relaxation time (ts), 
which could be an indication of enhanced spin-flip scattering of the electron spins in graphene in the 
presence of the molecular magnetic moments. The effect of the molecular layer is masked for low 
quality samples (low mobility), possibly due to dominance of Elliot-Yafet (EY) type spin relaxation 
mechanisms. 
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Graphene, one atom thick layer of sp^ carbon atoms, 
has potential for spintronic applications due to theoret¬ 
ically predicted high spin relaxation time (r^ 100 ns) 

and long spin diffusion length (A^ ~ 100 /im)^’^. These 
exceptional properties are attributed to negligible spin 
orbit coupling and weak hyperfine interaction due to the 
low atomic mass of carbon^. However, the maximum re¬ 
ported experimental values demonstrate of about 12 
/im^ for encapsulated graphene and about 2.7 ns for 
the hydrogenated graphene^, which although remarkable 
when compared with other metals and semiconductors, 
are still lower by more than an order in magnitude than 
the theoretically predicted values. A mismatch between 
theory and experiments suggests towards external fac¬ 
tors such as impurities/defects present near the graphene 
lattice, which dominate the spin relaxation process and 
result in a lower value for Ag. 

In order to probe the role of impurities on spin 
transport, one can systematically introduce them to 
graphene. In recent years, different research groups 
have demonstrated several ways of introducing impuri¬ 
ties (magnetic and non-magnetic) in graphene such as 
doping with adatoms, introducing defects and chemical- 
functionalization^”^^, each method introducing a differ¬ 
ent spin relaxation source. For example heavy metal 
atoms such as An can change the spin transport proper¬ 
ties in graphene via spin orbit coupling^On the other 
hand, light metal (Mg) ions can introduce charge impu¬ 
rity scattering of spins in graphene^^, although the ex¬ 
perimental study rules out the role of this mechanism^^. 
A significant change in the spin transport properties 
of graphene was reported in the presence of magnetic 
moments^^, which can be introduced via hydrogena¬ 


tion or by introducing vacancies in the graphene lattice. 
Remarkably, recent weak localization measurements on 
graphene^^ also show that magnetic impurities could play 
a key role in limiting the spin relaxation time in graphene. 
As proposed by Fabian et al.^^, if the localized moments 
are present at adatoms, they can act as spin hot spots and 
enhance the spin relaxation process via resonant scatter¬ 
ing. Therefore, the recent findings serve as an imperative 
to introduce magnetic impurities in graphene and inves¬ 
tigate their effect on the spin transport. 

Introducing the impurities via the methods described 
above may damage the graphene lattice and modify 
its electronic band structure^Alternatively, the self- 
assembly of molecular layers on graphene is a non¬ 
destructive way to functionalize the graphene surface 
and one can still tune the electronic properties of this 
two dimensional materiad^. Recently, Zhang et al. have 
reported the self-assembly of porphyrin ligand bound 
cobalt atoms (GoPP) on top of a graphene surface^^. 
Porphyrins are attached to graphene via weak Van der 
Waals interactions, while the cobalt atoms do not form 
any chemical bond with graphene in contrast to the direct 
deposition of metal atoms or ions as discussed above^^. 
Therefore, the self-assembly is not supposed to change 
the electronic properties of graphene significantly. On the 
other hand, cobalt atoms have an unpaired spin (S=l/2), 
which can act as a localized magnetic moment. 

In this work, we study spin transport properties of a 
CoPP-graphene system as a function of temperature, us¬ 
ing non-local spin-valve and Hanle spin precession mea¬ 
surements. After the self assembly of magnetic molecules, 
a reduced up to 50% with a lowered spin diffusion co¬ 
efficient Ds is obtained compared to the values for the 
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sample without functionalization (pristine sample). A 
pronounced effect of the molecular layer was observed for 
samples with high mobility and high diffusion coefficient, 
alluding to the sample quality playing an important role 
in determining the spin transport properties in graphene 
in contrast to previous studies^^. 

Graphene spin-valve devices are prepared using highly 
oriented pyrolytic graphite (HOPG), which has a very 
low amount of impurities (ZYA grade, supplier: SPI). 
Graphene is mechanically exfoliated on to a pre-cleaned 
Si/Si02 substrate (300 nm thick Si02), where n++ doped 
Si is used as a back gate electrode. Ferromagnetic (FM) 
contacts are patterned via electron beam lithography on 
the PMMA coated graphene flake. Then 0.4 nm of ti¬ 
tanium (Ti) is deposited in two steps, each step fol¬ 
lowed by oxidation to define a tunnel barrier, which is 
to overcome the conductivity mismatch problem^^. On 
top of the oxide barrier we deposit 35 nm of cobalt for 
the spin selective contacts. To prevent oxidation of the 
ferromagnetic electrodes, the contacts are covered with 
3 nm thick aluminum layer followed by the lift-off pro¬ 
cess. A lock-in amplifier detection technique is used to 
measure the charge and the spin transport properties of 
our samples. All the measurements are carried out us¬ 
ing a cryostat in vacuum (^ 1 x 10“^ mbar) at different 
temperatures between 4K and 300K. First, the sample 
is characterized in its pristine state. Afterwards, the 
magnetic impurities are added to the sample and the 
change in the charge and spin transport properties is 
measured. In order to equip graphene with magnetic 
molecules, a cobalt-porphyrin solution (cone. 0.56 mg/ml 
in tetradecane) is drop cast on top of the device and 
left to dry for 10 minutes. The residual porphyrin lay¬ 
ers on top are removed by rinsing the device with hexane 
(Fig. 1(b)). Since the exfoliated samples on the insulating 
Si 02 substrate are not big enough to perform scanning 
tunneling microscopy (STM), in order to confirm the self- 
assembly of porphyrins on graphene we perform STM on 
the large area GVD graphene-GoPP system. An STM 
image (Fig. 1(d)) of a GVD graphene sample (Si/Si02 
substrate) with the GoPP molecules on top confirms the 
self-assembly of cobalt-porphyrin molecules on graphene. 

We report the measurements for three samples, pre¬ 
pared under identical conditions. For discussion, they 
are labeled as sample A, B and G. A scanning electron 
microscope (SEM) image of sample A is shown in Fig¬ 
ure. 1(c). For the charge transport measurements, an 
alternating current (ac) is applied between contacts 1 
and 4 and the voltage is measured between contacts 2 
and 3 (Fig. 1(c)). In order to measure the carrier density 
dependence of the graphene resistivity (Dirac curve), we 
sweep the back-gate voltage. After the self-assembly of 
the GoPP molecules on the sample, the gate dependence 
is found to have a positive hysteresis at room temperature 
(inset Fig. 2), which alludes to a charge transfer process 
between graphene and the GoPP molecules^^’^^. At low 
temperatures charge states are frozen in the molecules 
and no hysteresis is observed. The field effect electron 
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FIG. 1. (a)Molecular structure of a cobalt bound porphyrin 

(CoPP) complex. Co^^ (in red circle) is the central atom in 
the complex, surrounded by the Porphyrin ligand. In the por¬ 
phyrin ring ‘-R’ represents a long chain alkyl group (-C 10 H 21 ), 
which is responsible for making weak Van der Waals interac¬ 
tion with graphene during the self-assembly, (b) Non-local 
measurement scheme for a graphene spin-valve. Graphene 
(in gray) with a self-assembly of cobalt-porphyrin molecules 
on top (cobalt magnetic moments in red), is probed with fer¬ 
romagnetic tunnel contacts (in blue), (c) Scanning electron 
micro-graph (SEM) of sample A. The distance between con¬ 
tacts 2-3 (transport channel) is 5 /xm. Outer contacts are 
chosen far enough from the inner ones, in order to make sure 
that they do not affect the spin transport, (d) A scanning 
tunneling microscopy (STM) image of CVD graphene func¬ 
tionalized with cobalt-porphyrin molecules on top (scan area 
39 nm^) on Si/Si02 substrate, which demonstrates an ordered 
self-assembly of the CoPP molecules on graphene. A bright 
spot in the image corresponds to the core of the porphyrin 
molecule. 


mobility jie for the pristine device is 7100 cm^V“^s“^, 
and for the GoPP device /Xg ^ 5000 cm^V“^s“^, both 
mobilities calculated at room temperature (RT) for a car¬ 
rier density ^ 10^^ cm“^. Gontact resistances (Rg) for 
all the samples were high enough (>1.5 kO) to be in the 
non-invasive regime as described in ref. 22. 

Eor the spin transport measurements, a four probe 
non-local detection scheme is used (Eig. 1(b)). This 
method allows us to decouple the charge and spin cur¬ 
rent paths and thus minimize the charge contribution 
to the detected spin signal {Rnl=^nl^ The spin- 
valve measurement is performed by first setting a high 
magnetic field ^ along the -y direction (Eig. 1(b)), so all 
the EM electrodes are magnetized along the field (parallel 
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FIG. 2. Resistivity as a function of gate voltage for the CoPP 
device (sample A) at different temperatures. Solid (dashed) 
lines in the plot correspond to forward (backward) sweeping 
direction of the back-gate voltage. The CoPP device shows 
hysteresis at room temperature (black curve), which disap¬ 
pears at low temperatures (blue curve). Hysteresis at RT in¬ 
dicates towards a charge transfer process between the CoPP 
molecules and graphene, which disappears at low tempera¬ 
tures due to freezing of the charged states in the molecules^^. 
A comparison between the Dirac measurement for the pristine 
and the CoPP state of the sample A is shown in the inset (at 
4 K). After functionalization, the sheet resistance increases 
near the charge neutrality point, which is not significant at 
high carrier densities. 


configuration). Then sweeping the field in the opposite 
direction, the electrodes reverse their magnetization at 
different fields depending on their coercivity, leading to 
an anti-parallel configuration between the inner injector 
and the detector electrodes, which appears as a switch 
in the non-local signal. At high magnetic field, all the 
electrodes are again magnetized in the same direction in 
the parallel configuration. The difference between the 
parallel and the anti-parallel signals is the spin-valve sig¬ 
nal ARnl- The outer contacts are chosen far away from 
the inner electrodes. In this way their influence on the 
measured spin signal is eliminated and we see only two 
distinct switches that correspond to the magnetization 
reversal of the inner injector and the inner detector. 

Spin-valve measurements for sample A before and af¬ 
ter the functionalization are shown in Fig. 3 at different 
temperatures. For both pristine and functionalized states 
of the sample, the spin-valve signal shows the switches 
corresponding to the contacts magnetization. However, 
after the functionalization, the signal magnitude is sig¬ 
nificantly reduced. At low temperature, the signal mag¬ 
nitude is increased for both the pristine and the CoPP 
devices (Fig. 3). 

In order to understand the effect of localized magnetic 
moments on spin transport in graphene, we refer to the 
exchange field model, explained by McCreary et al}^. In 
this model, an electron spin in graphene can interact with 


FIG. 3. Spin-valve measurements for sample A are shown in 
the positive x-axis for the pristine state and for the device af¬ 
ter the functionalization are shown in the negative x-axis. A 
strongly reduced spin-valve signal is observed after the func¬ 
tionalization. 


the magnetic moments via an exchange field which 
is the average exchange field produced by the localized 
moments. varies spatially and in time in a random 
fashion and influences locally the Larmor precession of 
the diffusing spins. The effect of varying precession fre¬ 
quencies at different locations resembles the D’yakonov- 
Perel mechanism of spin relaxation^^ and is responsible 
for an additional spin dephasing. In a spin-valve mea¬ 
surement, an enhanced relaxation (a reduced signal) is 
expected when the moments are randomized. As one 
starts applying an in-plane magnetic field, the magnetic 
moments try to align themselves along the field and their 
dephasing effect gets suppressed. This feature would ap¬ 
pear as a dip in the spin-valve signal. Within this picture, 
the spin relaxation rate by the fluctuating exchange field 
causing the dip is given by the following equation: 


1 _ 1 


( 1 ) 


where AB is the exchange field fluctuation magnitude, 
^e=2 is the gyromagnetic factor of the free electrons, jj^b 
is the Bohr magneton, h is the reduced Plank constant 
and Tc is the fluctuation correlation time^^. According to 
the formula above, the maximum relaxation (dip) in the 
spin-valve measurement should appear when Bapp=-Bex- 
Therefore the magnetic ordering of the localized moments 
affects the observation of the dip. For paramagnetic or¬ 
dering one would observe the dip around Bapp=0. On the 
other hand, for ferromagnetic ordering, there is a non¬ 
zero exchange field B^^ present {B^x 7 ^ 0 ) even when no 
external field is applied {Bapp=0). Now the dip would 
occur at finite external applied field and would exhibit 
hysteresis. 

For the measured spin-valve signal for the CoPP de¬ 
vice, we do not observe any dip, either around zero or 
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non-zero applied field. The only clear effect of intro¬ 
ducing the CoPP molecules is the reduced magnitude of 
the spin-valve signal. The observed behavior can be ex¬ 
plained by considering the magnetic moments playing the 
role of spin-flip scatterers in the transport channel, which 
enhance the spin relaxation process but do not produce 
a measurable effective exchange field. In order to con¬ 
firm if the source of the reduced spin signal is due to an 
enhanced spin relaxation rate, we now need to measure 
the spin transport parameters via Hanle spin precession 
measurements. 

Hanle precession measurement is a reliable tool to 
study the spin transport properties. Here, a magnetic 
field is applied perpendicular to the direction of the in¬ 
jected spins, which precess around this field ^ with Lar- 
mor frequency uJl = While sweeping the mag¬ 

netic field, due to the precession, spins can be reoriented 
to a direction opposite to the injected one, leading to 
a sign reversal in the spin signal. Simultaneously, they 
also dephase and result in a lower spin accumulation at 
higher fields. The Hanle precession can be fitted with the 
equation^^: 

^oo 1 t 

ARnloc / — e^^st cos{ujLt)e^dt^ ( 2 ) 

Jo {^irDst) 

where Dg is the spin diffusion coefficient, the spin re¬ 
laxation time and is the Larmor frequency. 

Referring back to the exchange field model, a Hanle 
measurement in the presence of an exchange magnetic 
field Bex by the magnetic moments would represent a 
spin precession due to a net field BappABex- The pre¬ 
cession can result in a narrower Hanle shape due to an 
enhanced g factor^d4,27 ^ paramagnetic ordering of 
the localized moments. Whereas for the case of ferromag¬ 
netic ordering, we would expect a shifted Hanle peak. 

Hanle precession curves for both pristine and CoPP 
devices are shown in Fig. 4. Here we show the purely 
spin dependent signal, obtained by subtracting the anti¬ 
parallel signal from the parallel signal and the result is 
fitted via Eq. 2. The plots have been normalized to the 
value at Bapp=0 for clearly demonstrating the change in 
the Hanle line shape. We observe two general trends for 
all measured samples. First, the Hanle curve becomes 
broader after the CoPP self-assembly. This is in con¬ 
trast to the expected narrowing of the Hanle curve in 
the presence of a paramagnetic exchange field according 
to the model described above. The observed broadening 
indicates a reduction of the spin relaxation time, in ac¬ 
cordance with our interpretation of the signal reduction 
in spin-valve measurements. Second, upon decreasing the 
temperature from RT down to 4 K we do not observe any 
significant narrowing of the Hanle line shape which could 
be interpreted as an enhanced g factor. On the contrary, 
the typical line widths and extracted spin lifetimes are 
not strongly dependent on temperature. 

A summary of the extracted spin parameters for all 
samples studied in this work is presented in Fig. 5. For 



B^(T) 



FIG. 4. a) Hanle measurements ((Rp-Rap)/2) for the pristine 
(black squares) and the CoPP state (red circles) at 4K (sample 
B). The corresponding fittings are plotted in line. The curves 
are normalized with respect to the signal at B=0. After the 
functionalization Hanle line shape is broadened, indicating a 
reduced spin relaxation time (ts). (b) Hanle measurements 
for sample B after self-assembly at RT and 4K. The curves 
are normalized. Broadening of the black curve (with square 
symbol) (RT) is dominated by the enhanced Dg. The spin 
relaxation time (r^) only changes from 100 ps (RT) to 112 ps 
(4K). 


sample A we observed the strongest effect of the molec¬ 
ular layer on the spin parameters. In its pristine state, 
the extracted spin relaxation time is in the range of 
300-400 ps for all the measured temperatures, with a cor¬ 
responding spin diffusion length Xs{= \/DsTs) 3-4.5 /im. 
On the other hand, after self-assembly sample A exhib¬ 
ited a strongly reduced Tq in the range 100-200 ps and 
a correspondingly lower Ag ^ 2-2.5 /im. Interestingly, 
we did not observe any significant temperature depen¬ 
dence for the extracted in the measured temperature 
range, which would have otherwise been expected due 
to the presence of an effective exchange field via local¬ 
ized molecular paramagnetic moments^^’^^. Therefore 
the added magnetic molecules seem to only increase the 
spin relaxation rate via the introduction of more spin-flip 
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FIG. 5. A summary of (a) Ts and (b) Ds, extracted from 
Hanle analysis, for sample A (square+line), B (triangle) and 
C (circle) before (black) and after (red) the functionalization. 
Black data corresponds to the pristine and red data is for the 
CoPP state of the samples. A reduced Ts and Dg were ob¬ 
served for the samples after the functionalization with a weak 
temperature dependence, which rules out any exchange cou¬ 
pling between the localized magnetic moments and the elec¬ 
tron spins in graphene^ ^ and indicates towards an enhanced 
spin-flip process, where the present magnetic moments play 
only the role of spin-flip scatterers. The effect of the molec¬ 
ular layer is determined by the sample quality (fie, D) in the 
pristine state. Since sample A and C have higher mobility and 
diffusion coefficient, Tg is highly reduced for these samples af¬ 
ter the functionalization. Sample B, having lower mobility 
did not show any significant change in Tg. 


scattering events. 

Furthermore, we also observed a minor reduction of the 
extracted spin diffusion coefficient Dg after self-assembly, 
consistent with the observed reduction in mobility as dis¬ 
cussed in Fig. 2. Note that the reliability of a Hanle fit 
is typically established by comparing the agreement be¬ 
tween the extracted spin diffusion coefficient Dg with the 
charge diffusion coefficient The latter can be 

independently calculated via the resistivity of the sample 
at a known density of states u using the Einstein relation 
Dc = XjRgq^v. In the absence of electron-electron inter- 


4K 

RT 


pristine 

CoPP 

pristine 

CoPP 

Dc Ds 

Dc Ds 

Dc Ds 

Dc Ds 

A 

0.052 0.048 

0.039 0.034 

0.100 0.037 

0.050 0.027 

B 

0.010 0.010 

0.011 0.010 

0.014 0.012 

0.020 0.020 

C 

- - 

0.12 

0.085 

0.063 


TABLE I. A summary of and Dg (units in vni^ js) for sam¬ 
ple A, B and C, before (pristine) and after (CoPP) function¬ 
alization. For all the samples, Dc and Dg are approximately 
in the similar order. For sample A, Dc in the pristine state 
is found around 0.05 vr? js. We also sometimes observed an 
asymmetry in the Dirac curve at different temperatures. This 
asymmetry rises due to contact induced doping at different 
regions^®, resulting in a different value for Dc at different 
temperatures. 


action D^ and Dg should match^^. As shown in Table. I 
both parameters are in a reasonable agreement, confirm¬ 
ing the validity of the the Hanle analysis. 

It is worth mentioning that the earlier work of molecu¬ 
lar doping on graphene^ ^ did not exhibit any measurable 
change in the spin transport properties of graphene, while 
the charge transport properties were modified. However, 
we find that sample quality, as determined by the mag¬ 
nitude of the diffusion coefficient [D^^ Dg) or electronic 
mobility (/ie), plays an important role on the influence 
that the cobalt-porphyrin molecular layer exerts on the 
spin transport parameters. For example, for sample B 
we do not observe a significant change in Dg and Tg after 
self-assembly. This reduced sensitivity can be attributed 
to its low mobility (and diffusion coefficient) which in the 
pristine state was ^2000 cm^V“^s“^, almost a factor of 
3.5 times lower than for sample A. On the other hand, 
for sample C which had a comparatively better quality 
{Dg ^0.1 m^/s) we again observed a significant reduc¬ 
tion of 30 % on the spin relaxation time, confirming our 
initial observations. 

A significant reduction in r^, with a simultaneous mod¬ 
erate reduction in Dg, is inconsistent with the picture 
of localized magnetic moments creating an effective ex¬ 
change field as discussed above^^’^^, or a model where 
localized states act as spin reservoirs^^ . Both models 
imply a significant increase of the extracted Tg and a 
proportionally reduced Dg, which can be understood via 
an enhanced g factor and the symmetry of the Hanle 
equation^’^^. Furthermore, both models are also ex¬ 
pected to show a strong temperature dependence, which 
is not observed here. 

The reduction on the spin transport parameters indi¬ 
cates that the main role of the Co-porphyrin molecular 
layer is to act as an extra source of spin-flip scattering. 
This interpretation is consistent with the lack of sensitiv¬ 
ity to the molecular layer by low-quality samples, where 
the initial spin relaxation rate was already large and 
therefore masks the relaxation process introduced by the 
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molecular layer. In addition, the concomitant reduction 
in Dg and observed can be partially understood by the 
enhanced momentum scattering introduced by the molec¬ 
ular layer, since in single layer graphene the leading spin 
relaxation mechanism is of the Elliot-Yafet type, which 
results in the proportionality relation Tg oc This 

observation is interesting, since previous experiments rule 
out the role of mobility dependence of or seem to ob¬ 
serve an opposite relation between Tg and fie i.e. higher 
spin lifetime for lower mobility samples^^. 

To summarize, we observe a change both in the charge 
and spin transport properties of graphene in the presence 
of cobalt-porphyrin molecules. In the charge transport 
measurements, we observe an increase in the graphene 
sheet resistance after functionalization due to their in¬ 
teraction with graphene via weak Van der Waals forces. 
For the spin transport measurements we observe lower 
values of Tg and for the CoPP-graphene system com¬ 
pared to the pristine one. The measurements are not 


strongly temperature dependent suggesting that the ad¬ 
ditional spin relaxation can be due to the enhanced spin 
scattering from the magnetic moments, accompanied by 
the other relaxation sources. The changes are also sen¬ 
sitive to the sample quality {D^fie) in the pristine state 
and are masked for a lower value of the mobility or dif¬ 
fusion coefficient, indicating also the presence of Elliot- 
Yafet type spin relaxation mechanism. 
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